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A Convex Approach to Hydrodynamic Analysis 

Mohamadreza Ahmadi, Giorgio Valmorbida, Antonis Papachristodoulou 


Abstract —We study stability and input-state analysis of three 
dimensional (3D) incompressible, viscous flows with invariance 
in one direction. By taking advantage of this invariance prop¬ 
erty, we propose a class of Lyapunov and storage functionals. 
We then consider exponential stability, induced £ 2 -norms, and 
input-to-state stability (ISS). For streamwise constant flows, 
we formulate conditions based on matrix inequalities. We 
show that in the case of polynomial laminar flow profiles the 
matrix inequalities can be checked via convex optimization. 
The proposed method is illustrated by an example of rotating 
Couette flow. 


I. INTRODUCTION 

The dynamics of incompressible fluid flows is described 
by a set of nonlinear partial differential equations known as 
the Navier-Stokes equations. The properties of such flows 
are then characerized in terms of a dimensionless parameter 
Re called the Reynolds number. Experiments show that 
many flows have a critical Reynolds number Rec below 
which global stability is ensured. However, spectrum analysis 
of the linearized Navier-Stokes equations, considering only 
infinitesimal perturbations, predicts a linear stability limit 
Rer which upper-bounds Rec [1]. On the other hand, the 
bounds using energy methods ReE, the limiting value for 
which the energy of arbitrary large perturbations decreases 
monotonically, are much below Rec [2], For example, 
ReE = 32.6 [3], ReL = oo [4] and Rec = 350 [5] for 
3D Couette flow. 

The discrepancy between ReL and Rec have long been 
attributed to the eigenvalues analysis approach [6], citing a 
phenomenon called transient growth ; i.e., although the per¬ 
turbations to the linearized Navier-Stokes equation are stable, 
they undergo high amplitude transient amplifications that 
steer the trajectories out of the region of linearization. This 
has led to studying the resolvent operator or e-pseudospectra 
based on the general solution to the linearized Navier- 
Stokes equations [7]. Another method for studying stability is 
based on spectral truncation of the Navier-Stokes equations 
into an ODE system. Recently in [8], [9], a method was 
proposed based on keeping a number of modes from Galerkin 
expansion and bounding the energy of the remaining modes. 
However, these bounds on Rec turn out to be conservative. 

Since the seminal paper by Reynolds [10], it was observed 
that external excitations and body forces play an important 
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role in flow instabilities. Mechanisms such as energy ampli¬ 
fication of external excitations have shown to be crucial in 
understanding transition to turbulence [2]. Energy amplifica¬ 
tion of stochastic forcings to the linearized Navier-Stokes 
equations in parallel channel flows was studied in [11], 
[12], In [12], using the linearized Navier-Stokes equation, 
it was shown analytically, through the calculation of traces 
of operator Lyapunov equations, that the 77, 2 -norm from 
streamwise constant excitations to perturbation velocities is 
proportional to Re 3 . The 0(Re 3 ) amplification mechanism 
of the linearized Navier-Stokes equation was verified in [13] 
and [14], where the influence of each component of the body 
forces was calculated in terms of R 2 and 77 00 -norms. Input- 
output analysis of a model of plane Couette flow was carried 
out in [15] to study the nonlinear mechanisms associated with 
turbulence. In another vein, an input-state analysis method 
for the linearized Navier-Stokes equation by calculating the 
spatio-temporal impulse responses was given in [16]. 


In this paper, we study the stability and input-state proper¬ 
ties of incompressible, viscous fluid flows. We study input- 
state properties such as induced £ 2 -norms from body forces 
to perturbation velocities and ISS. In particular, we consider 
flows with invariance in one of the three spatial coordinates. 
For such flows, we formulate a suitable structure as a 
Lyapunov/storage functional. Then, based on these function¬ 
als, for streamwise constant flows, we propose conditions 
based on matrix inequalities. In the case of polynomial 
laminar velocity profiles, e.g. Couette and Poiseuille flows, 
these inequalities can be checked via convex optimization 
using available computational tools. The proposed method is 
applied to the analysis problem of a rotating Couette flow. 


The paper is organized as follows. The next section 
presents some preliminary results. In Section m we 
formulate the Lyapunov/storage functional stmcture. 
Section m is concerned with the convex formulation 
for streamwise constant flows. The proposed method is 
illustrated by studying an example of a model of rotating 
Couette flow in Section [V] Finally, Section VI concludes 
the paper and provides directions for future research. 


Notation: The n-dimensional Euclidean space is denoted by 
R n . The nxn identity matrix is denoted by I nxn . A domain 
U C R” is a connected, open subset of R", and is the 
closure of set f l. The boundary Oil of set S2 is defined as 
fl\Sl with \ denoting set subtraction. The space ofp-th power 
integrable functions defined over U is denoted £[', endowed 
with the norm ||(-)||£p = (//q(‘) p dfl) p , for 1 < p < oo, 
and IKOlUff = suPzen KOI, for P = oo. Also, we denote 



by C^ to T j n , with to > 0, the space of square integrable 
functions in x £ Cl and t £ [to,T] with the norm 



The space of fc-times continuous differentiable functions 
defined on f l is denoted by C k (Cl). If p £ C 1 , then d Xl p 
is used to denote the derivative of p with respect to variable 
Xi, i.e. d Xl := -^T-. A continuous strictly increasing function 
k : [0,a) p —> K>o, satisfying k( 0) = 0, belongs to class 
1C. If a = oo and linx^oo k{x) = oo, k belongs to class 
/Coo. The unit vector in direction Xi is denoted by T 2 ,. For a 
scalar function v , X7v = d Xi v~ifi i denotes the gradient and 
X7 2 v = dfv denotes the Laplacian. For a vector valued 
function w = Y, w i e i, the divergence V ■ w is given by 
V • w = J2i d Xi Wi. 

II. Preliminaries 

A. Flow Model 

We consider incompressible, viscous flows with invariance 
in one of the directions^ x m , m £ {1,2,3}, i.e., d Xm = 0. 
Let I = {1, 2, 3} — {m\. The flow dynamics is described by 
the Navier-Stokes equations, given by 

d t u = — X/ 2 u — u ■ X/u — X/p + Fu + d, 

Re 

0 = V-m, (1) 

where t > 0, F £ M 3x3 , and x £ !1 = il t x Cfi clxl 
with x = (xi, Xj)', i,j £ I being the spatial coordinates. The 
dependent variable d(f,x) = [rii(f,x) d 2 (f,x) d 3 (f,x)] 

is the input vector representing exogenous excitations or 
body forces, tt(/,x) = [iti(/,x) U 2 (t,x) U 3 (f,x)] is the 

velocity vector, and p(t, x) is the pressure. 

We consider perturbations (it, p) to the stationary flow 
(U,P). That is, 

u = u + U, p = p + P, (2) 

where (U,P) satisfy 

0 = —X/ 2 U - U ■ X/U - VP + FU, 

Re 

0 = V • U. (3) 

Substituting ([2]i in 0 and using ([3]), we obtain the pertur¬ 
bation dynamics 

dtu = —V 2 tt — u ■ X7u — U ■ X7u — u ■ X7U 

Re 

—X/p + Fu + d, 

0 = V-u. (4) 

In this paper, we concentrate on perturbations with no¬ 
slip boundary conditions it|gf) = 0 and periodic boundary 
conditions. 

1 Invariance in one direction is a common assumption in the case of several 
fluid models, namely, Couette flow, Poiseuille flow, Taylor-Couette flow, etc. 


B. Stability and Input-to-State Analysis 

In this section, we briefly review a number of definitions 
and results from [17] and [18]. 

Definition 1 (Exponential Stability): The stationary solu¬ 
tion (0,po) of 0 with d = 0 is exponentially stable in tZ 2 ,, 
if there exists a A > 0, such that for all i > 0 

|Kf,x)|| 2 a <M0,x)||^e— At . (5) 

That is, system 0 converges to the laminar flow (U, P) as 
in 0. 

Definition 2 (input-to-State Properties): 

A. Induced C 2 -norm Boundedness: For some 

rji > 0 , * = 1 , 2 , 3 , 

3 

X )H £ fo,=o),n < Z^ i ll di (*’ X )ll£fo,oo ) ,n (6) 

i—1 

subject to zero initial conditions u{ 0,x) = 0, Vx £ Cl. 

B. Input-to-State Stability: For some scalar y: > 0, func¬ 
tions ft, /?, x £ /Coo. an d cr £ /C, it holds that 

ll«(f> x )ll/;2 < ^{er^x (lM0,x)|| £ 2j j 

+ p( sup ([[ <7(|d(r,x)|) dfi)Y (7) 
\re[o,t) JJn / 

for all t > 0. 

Remark 1: Due to nonlinear dynamics, the actual induced 
£ 2 -norms of system 0 are nonlinear functions of ||d||£ 2 . 
The quantities rji, i = 1,2,3 provide upper-bounds on the 
actual induced £ 2 -norms. 

Remark 2: The ISS property 0 implies the exponential 
convergence to the laminar flow (U, P) in C 2 , when d = 0. 
Moreover, as t —>■ oo, we obtain 

^ ll M (f> x )ll £2 < ^ fj ^ ||cr(|d(/,x)|)|| £ - oo) d Clj 

</3 (^0'(l|d(/, x )|U~ oo) ) dfl) , (8) 

wherein, the fact that cr, ft £ 1C is used. Hence, as long as the 
external excitations or body forces d are bounded in 
(this encompasses persistent excitations), the perturbation 
velocities u are bounded in Cq sense. 

The next result converts the tests for exponential stability, 
induced £ 2 -norm boundedness, and ISS into the existence 
problem of a Lyapunov or a storage functional satisfying a 
set of inequalities. 

Theorem 1: Consider perturbation model {4}. If there exist 
a positive definite Lyapunov functional V(u) and a pos¬ 
itive semidefinite storage functional S(u). positive scalars 

{»?*}ie{i,2,3}. {ci}te{i,2,3}> t/t* and functions /3 i,/ 3 2 £ /Coo, 

cr £ 1C, such that 
I) when d = 0, 

Ci|Ml £2 < V(u) < c 2 ||u||£ 2 , (9) 

d t V{u) < -c 3 ||tt||£ 2 , (10) 


II) 


d t S(u) < 


u'u dfl 



vl oo 
o vl o 
0 0 1)1 


ddfl , ( 11 ) 


III) 


PidWlcf ,) < S{u) < 02{\\u\\ C 2 n ), 


( 12 ) 


d t S(u) < —ipS(u) A 



cr(|d(t, x )|) dU, 


(13) 


for all t > 0 , then, respectively, system 0 

I) is exponentially stable, 

II) has induced C 2 -norm upper-bounds rj, , i = 1,2,3 as 
in <j 6 }, 

III) is ISS and satisfies (|7} with \ = P 2 , P(') = Pf 1 0 2(-) 
and f}{-) =/3f 1 o |(-). 

Proof: This is a direct application of Theorem 1 in [17] 
and Theorem 1 in [18]. ■ 


III. Lyapunov and Storage Functionals 
for Fluid Flows 

In this section, we derive classes of Lyapunov and storage 
functionals suitable for analysis of system 0 subject to 
invariance in one of the three spatial coordinates. In the 
following, we adopt Einstein’s multi-index notation over 
index j, that is the sum over repeated indices j, e.g., 

Vjd XjUj = V Vjd^Uj. 

The perturbation model {4} can be re-written as 

Of u.; V u, Ujd x .Ui Ujd x ,Ui 

Re 3 3 

m j 0 . t , u„ 0 T ; 7 p A Fijiij T d !: , 

0 = d X] u r (14) 

where i,j £ {1,2,3} and F.\j is the (i, j) entry of F. 

The next theorem states, under which Lyapunov/storage 
functional structure, the time derivative of the Lya¬ 
punov/storage functional takes the form of a quadratic form 
in dependent variables u and their spatial derivatives, by 
removing the nonlinear convection and pressure terms. 

Proposition 1: Consider the perturbation model ( [14} sub¬ 
ject to periodic or no-slip boundary conditions u\gn = 0 . 
Assume ( [14} is invariant with respect to x m , m £ {1, 2, 3}. 
Let I = {1, 2, 3} — {in} and 


V{u) 



k m 0 0 
0 kt 0 
0 0 kj 


u dSl 


3 


kiUi(t,x) 2 d£7, (15) 


where k, = kj for i,j £ I, be a candidate Lyapunov or 
storage functional. Then, the time derivative of C 0 satisfies 




i—1 


(~R^ U ^ + U i Uidx i Ui 


+ UjUid x .Ui - UiFijUj ) dS7, (16) 


Proof: The time derivative of Lyapunov functional ( [15} 
along the solutions of ([14} can be computed as 


d t V(u) = ^ JJ ki ( - UjU z d Xj i 

~UjUid x Ui — UjUjd x Ui — Uid Xi p+UiFijUj+Ujdi I df 2 . 


(17) 


Consider Jf n kiUjUid Xj Ui dfl. Using the boundary condi¬ 
tions, integration by parts and the incompressibility condition 
d X: . Uj = 0 , we obtain 



k, Ujllflr, 11 , dfl = 


1 

2 



j k,u :) v.,j | d Xi 


Qi 

kiU 2 (d Xj uj) dfl = 0 . (18) 


At this point, consider the pressure terms ff n kiUid Xi p dfl. 
Without loss of generality, we consider invariance in * 1 , 
which yields 


JJ^ {k 2 u 2 d X2 p A k 3 u 3 d Xa p) dU 

= / (k 2 u 2 p)\gn 2 dx 3 + / {k 3 u 3 p)\gn 3 dx 2 
J Q 3 J Q 2 

// {k 2 d X2 u 2 p + k 3 d X3 u 3 p) dU 

J Jfl 

= ~ f [ ( k 2 d X2 u 2 A k 3 d Xa u 3 )p dU, (19) 
J J n 


where in the first equality above integration by parts and 
in the second inequality the boundary conditions are used. 
Then, if k 2 = k 3 , using the incompressibility condition 
T d X3 U 3 0, {±9} equals zero. Therefore, the time 
derivative of the Lyapunov/storage functional ([17} is modi¬ 
fied to 


()i L (/h) — } 7 JJ ^1 (/, ll 7 () :r ^ IL, 


- UjUid x ,Ui A UiF^Uj A Uidi d£7. (20) 


Integrating by parts the UiX/ 2 Ui term and using the boundary 
conditions, we get 


d t V(u) = ^ JJ ki Jjf( dx i u i) 2 ~ UjUid Xj Ui 

UjUiO X j Ui A UiFijUj A UidJ d- 1. (21) 

Applying Poincare inequality (Lemma 1 in Appendix A) 
to ( | 2 l} , we obtain ©■ ■ 

Remark 3: A special case of GD was used in [19] to 
study the stability of viscous fluid flows. 

Remark 4: In the sequel, we use structure {B} as a 
Lyapunov functional when studying stability and as a storage 
functional when studying input-to-state properties. 


where C > 0. 







Remark 5: There are several estimates for the optimal 
Poincare constant. The optimal constant we use in this paper 
is 

C(n > = mr <22) 

where D(fl) is the diameter of the domain Cl [20], 

The next corollary proposes conditions under which prop¬ 
erties such as stability, input-state induced C 2 bounds and 
ISS can be inferred for the flow described by ( [T4| . 

Corollary 1: Consider the flow described by ( fl4| ) subject 
to periodic or no-slip boundary conditions u | ,-)o = 0. As¬ 
sume the flow is invariant with respect to x m , to £ {1, 2, 3}. 
Let I = {1,2,3} — {to}. If there exist positive constants 
hi, i = 1,2,3, with ki = kj , i,j £ /, positive scalars 
{A}ie{ i. 2 , 3 }, Mie{i, 2 , 3 }, and a £ 1C such that 

I) when d = 0, 

+UjUid Xj Ui — UiFijUj^J dfi > 0 (23) 

II) 


mw 



T kjUjiijd X j Uj kiUiFijUj 

HI) 


+ kiUjUid Xj Ui 
kiUidi + ijld 2 ) 


<m > 0 

(24) 


T JJ ^ i~i^~ ~ ^ ik J + kiUjuAjUi 

+ k, uj iijd x I 'i - kiUiFijUj 

- huidi + cr(|di|, \d 2 \, dfl > 0 (25) 


Then, 

I) perturbation velocities given by © are exponentially 
stable. Therefore, the flow converges to the laminar flow 
exponentially. 

II) under zero perturbation initial conditions u (0.x) = 0, 
the induced £ 2 norm from inputs to perturbation velocities 
is bounded by i £ {1, 2, 3} as in ([ 6 }. 

III) the perturbation velocities described by ( [14} are ISS in 
the sense of 0- 

Proof: Each item is proven as follows. 

I) Considering Lyapunov functional ( p~5] l, inequality 0 is 
satisfied with c\ = min ie r 1 2j 3 } ki and c 2 = maxj e / 1)2 , 3 } ki. 
Re-arranging the terms in (|23| gives 




(cm 2 


\ Re 


Uj U j'Hidxj'Ui T U jUi d Xj Ui 


—UiFijUj — Uidi ) dfl < 0 . (26) 


Then, from Proposition [T] we infer that, for d = 0, 
dtV(u) < 0. By continuity, we infer that there exists C 3 > 0 
such that © holds. Then, form Item I in Theorem 1, we 
infer that the perturbation velocities are exponentially stable. 

II) Re-arranging terms in (|24} yields 


^ ' ki J J ^ ^ Hi T U jUjd X jUj T VjjVjidx-Ui 


— UiFijUj — Uidi ) dS 2 


i—1 


< E // ^ dfi +E // dn (27 > 


2 — 1 


Then, from © in Proposition [T] we deduce that 


d t V(u) <-J2 u2 i cm + y dS1 - 


From Item II in Theorem 1, we infer that, under zero initial 
conditions, the perturbation velocities satisfy 0 

III) Adopting © as a storage functional, © is sat¬ 
isfied with /?!(•) = min l 6 {li 2 i 3 } fc ^-) 2 and /3 2 (-) = 
maxjg^ 2 3 } ki(-) 2 . Re-arranging the terms in ( |25) , we ob¬ 
tain 


kiCm) u 2 + k u q + kiUjUid x Ui 
Re 3 3 

— kiUiFijUj — kiUidi ) dll 


E 

2=1 


<~y^i // hu 2 dfl+ // cr(|di|, \d 2 \, |d 3 |) dCl 

JJ n JJn 

(28) 


From © in Proposition |T] it follows that 

d t V(u) < -tpV(u) + [[ a(\di\,\d 2 \,\d 3 \) dfl, (29) 
J Jn 

with ip = niiii ;f= { i 2 3 | tpi. Then, from Item III in Theorem 1, 
we infer that the perturbation velocities satisfy the ISS 
property 0- ■ 

IV. Convex Formulation for Streamwise 
Constant Flows 

To present a convex method for checking the conditions in 
Corollary [T] we restrict our attention to streamwise constant 
flows in x m -direction with laminar flow U = U m (x)lt m . 

Corollary 2: Consider the perturbation dynamics given 
by ( |T4| ). Assume streamwise invariance in .xy r , -direction with 
laminar flow U = U m (x)l( m where m £ {1,2,3}. Let I = 
{1, 2, 3} —{to}. If there exist positive constants {fc;}je{i, 2 . 3 } 
with k p = k q , p,q £ I, {m}ie { 1 , 2 , 3 }, 1 , 2 , 3 }, and 

functions {<tz}z 6 {i j2i 3 } such that 
I) © holds, 









> 0, i,j G I, j, xefi. 


M(x) = 


— F ) k 
V Re 1 771171 ) rx " m 
km (d Xj U m (x) - F mj ) -kj Fjm 

2 

km(O x . Um{x)-Fmi)-kiFim 


km(d x .Um(x)- Fmj ) - kj Fj r 


{Re Fjjjhj 


km(d Xi Um{*)-Fmi)-kiFim- 


k-jFn, 


(F_ _ L\. 
KRe r *» 


( 30 ) 


II) 


iV(x) = 





k m 

2 

0 

0 

M(x) — I 3x 3 

0 

kj 

2 

0 


0 

0 

ki 

2 

-¥ o o 

An 

0 

0 

O 

■irh 

i 

o 

0 

A 

0 

1 - 

o 

o 

1 

0 

0 

A 


for ^ j and xe!!, 

HI) oi (x) >0, x e fl, l e {1, 2 , 3} and 


P(x) = 




km 

2 

0 

0 

M(x) - Q 


0 

kj 

~~2 

0 



0 

0 

ki 

2 

-Min Q 

2 u 

0 

OVn(x) 

0 

0 

jP|« 

1 

o 

0 

0 

<7j( x ) 

0 

0 0 

ki 

2 

0 

0 

o'i(x) 


> 0 , 


1pm k m 0 0 

0 i/jjkj 0 
0 0 'ipiki 


for i, j £ I,i 7 ^ j and x £ Q, where Q = 

Then, it follows that 

I) the perturbation velocities are exponentially stable, 

II) subject to zero initial conditions, the induced C 2 norm 
from inputs to perturbation velocities is bounded by r/,, 
i = 1 , 2 ,3 as in ([ 6 |, 

III) the perturbation velocities are ISS in the sense of ([7]) 
with cr(|d|) = £i=i°i( x ) d i- 

Proof: The proof is straightforward and follows from 
computing conditions ((23), ((27), and ( |25| considering x m - 
invariance, the laminar flow U = U(dfe m , and u(|d|) = 
Si=i 0 i(x)d? ■ Since the flow is a; m -invariant and the lam¬ 
inar flow is given by U = Umltm, UjUjd x Ui = 0, 
* = 1,2,3. 


I) Inequality 

A = 


( cm 

V Re 

(cm 

^ Re 

(cm 


>0, (31) 


is given by 


Fa j kj '//, /i,( kjkjj)iij Ui(kjFi m ) 

n I kj u j — u j(kjFji)ui — Uj(kjFj m )u 


k rn e. vl 


I ^ hum j ^m^rn 


n L 


+ UmidxJJm Rmi)Ri 
+ u m {d Xj Um - Fmj) ) df 2 > 0 (33) 


(34) 


for i, j G /, i f j, which can be rewritten as 


M(x) 


dS 2 > 0 . 


with M(x) given in ( |3()1 ). Therefore, if ( [30) is satisfied, ( |34) 
also holds and from Item I in Corollary [I] we infer that the 
perturbation velocities are exponentially stable. 

II) Inequality (|27j) is changed to 


A + I / (kiUidi + kjUjdj + k m u m d m ) dfi 
In 


[u\ + A + An) dll 


+ / / (»h di + Vjdj + V m d m ) > 0, (35) 


for i,j £ I, i ^ j, which can be rewritten as 


an > 0 , 


( 32 ) 

IL 

r Um - 
Uj 

Ui 

dm 

dj 

/ 

iV(x) 

_ Um -i 

Uj 

Ui 

d m 

dj 



L di J 


L di -1 


(36) 


where N is defined in ( |3T) . Consequently, if a is satisfied 
for all x £ fl, ( [36) holds and from Item II in Corollary [I] 
we infer that, subject to zero initial conditions, the induced 
C 2 norm from inputs to perturbation velocities is bounded 
by rji, i = 1,2,3 as in ([ 6 ). 

Ill) The proof follows the same lines as the proof of Item II 
above. ■ 

In the case that U m (x) is a polynomial in x, inequali¬ 
ties a, a- and ( |32) are polynomial matrix inequalities 
that should be checked for all x £ Pi. If the set (2 is a 
semi-algebraic set then these inequalities can be cast as a 
sum-of-squares (SOS) program (see Appendix B) by using 
Putinar’s Positivstellensatz theorem [21, Theorem 2.14]. 

Remark 6: In order to find upper-bounds on the induced 
£ 2 -norm from the body forces (di, d 2 ,df) to the perturbation 
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Fig. 1: Schematic of the rotating Couette flow geometry. 



Fig. 2: Estimated critical Reynolds numbers Re in terms of 
Ro for rotating Couette flow. 


velocities u, we solve the following optimization problem 

minki^iVi+vZ + vl) 

subject to N(x) > 0, ki,kj > 0, i,j £ I. (37) 
In the next section, we consider the analysis of the rotating 
Couette flow, which illustrate the proposed results. 


V. Example: Rotating Couette Flow 


We consider the flow of viscous fluid between two co¬ 
axial cylinders, where the gap between the cylinders is 
much smaller than their radii. In this setting, the flow 
can be schematically illustrated as in Figure |T] The axis 
of rotation is parallel to ir 3 -axis and the circumferential 
direction corresponds to x.[ -axis. Then, the dynamics of 
the perturbation velocities is described by 0. The flow 
is assumed to be invariant with respect to X\ {d Xl = 0) 
and periodic in x 3 with period L. Therefore, O = 
{{x 2 ,xs) | (X 2 ,X 3 ) £ [—1,1] x [0,L]}. The laminar flow is 
given by U = ( X 2 ,0, 0)' = x 2 ~e s ’i and P = Pq. In addition. 


F = 


0 

—Ro 
0 


Ro 0 
0 0 
0 0 


where Ro £ [0,1] is a parameter representing the 

Coriolis forced We consider no-slip boundary conditions 

«IL =_1 =0 and u(t, x 2 , x 3 ) = u(t, x 2 , x 3 + L). The 

2 

Poincare constant is then given by C = L2+2 2 • 

Notice that the cases Ro = 0,1 correspond to the Couette 
flow. Thus, the obtained results for rotating Couette flow can 
be applied to the Couette flow in special cases, as well. We 
are interested in finding estimates of the critical Reynolds 
number Rec using the following Lyapunov functional 


V(u) = 


f L f 1 

ni ‘ 

/ 

o 

o 


' u-i ■ 


u 2 


0 fc 2 0 


U 2 

/ 0 J- 1 

U3 


0 0 k 2 


U3 


dx 2 d.T 3 , 


which is the same as Lyapunov functional ( [15] ) considering 
invariance with respect to x\. 


1 That is, Ro = 0 (Ro = 1) corresponds to the case where only the outer 
(inner) cylinder is rotating and Ro = 0.5 is the case where both cylinders 
are rotating with the same velocity but in opposite direction. 


For stability analysis, we need to check inequality ( |30l ) 
according to Item I in Corollary [2] For this flow (to = 1, j = 
2, i = 3), we have 


M = 


fci c 

Re / x 
k 2 Ro—k\ {Ro— 1) 

2 


o 


k 2 Ro—k\ {Ro— 1) 
2 

k 2 C 

Re 

o 


o 

0 

k 2 C 
Re - 


>0 (38) 


This is a linear matrix inequality (LMI) feasibility problem 
with decision variables k \, /c 2 > 0 . 

To find estimates of Rec in the case of Couette Flow 
Ro = 0, applying Schur complement theorem [22, p. 650] 
to ( [38] ), we have 


kiC 

Re 





which yields the incqualit\£] 




(39) 


This implies that the Couette flow is stable for all Re. Hence, 
for Couette flow, Rec = oo obtained using Lyapunov func¬ 
tional ( p~5] > coincides with linear stability limit Rcl = oo [4]. 

Let L = 7 T. Figure [2] illustrates the estimated critical 
Reynolds numbers Rec as a function of Ro obtained from 
solving the LMI ( |38| > and performing a line search over Re. 
Notice that for the cases Ro = 0,1 the flow is stable for all 
Reynolds numbers. 

For induced £ 2 -norm analysis, we apply inequality © 
which for this particular flow is given by the following LMI 
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M — I 3x3 
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£ 1 " 
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.£|<n 

r]f 0 0 
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0 0 rg 


3 For Ro = 1, we can similarly obtain -A- > 
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Fig. 3: Upper bounds on induced £ 2 -norms from d to perturbation velocities u of rotating Couette flow for different Reynolds 
numbers: Re = 2 (left). Re = 2.8 (middle), and Re = 2.83 (right). 


with M as in ( |38| . 

Figure [3] depicts the obtained results for three different 
Reynolds numbers. As the Reynolds number approaches the 
estimated Rec for Ro = 0.5, the upper-bounds on the 
induced £ 2 -norm from the body forces d to perturbation 
velocities u increases dramatically. 

The obtained upper-bounds on the induced £ 2 -norm for 
Couette flow Ro = 0, are also given in Figure [4] Since the 
flow is stable for all Reynolds numbers, the induced C 2 - 
norms keep increasing with Reynolds number. The obtained 
upper-bounds depicted in Figure [4] are consistent with Corol¬ 
lary 2 and Corollary 4 in [13] and Theorem 1 in [12], wherein 
it was demonstrated that r]f ex 0(Re), and oc 0(Re 3 ) 
for Couette flow. 

In order to check the ISS property, we check inequal¬ 
ity ( |32| i from Corollary [2] for the rotating Couette flow under 
study, i.e.. 



o 

o 

1 

M-Q 

o 

£1" 

1 

o 


0 0 - 'f 

O 

o 

1 

3 

o 

o 

o 

1 

o 

0 tJ2 0 

0 0 

0 0 (73 


ki'ifji 0 0 

0 & 2’02 0 

0 0 k 2 


We fix 


with M given in ( [38] ) and Q = 

ipi = 10 -4 , i = 1,2,3 and L = 2n. Figure |5J depicts the 
maximum Reynolds number for which ISS certificates could 
be found Reiss and stability critical Reynold’s numbers 
Rec as a function of Ro. It appears that for Ro £ (0,1) these 
two quantities coincide. However, for the case of Couette 
flow Ro = 0,1, we obtain Reiss = 316 and Rec = oo. 
The quantity Reiss = 316 is the closest estimate to the 
empirical Reynolds number Re ~ 350 [5]. 


VI. CONCLUSIONS AND FUTURE WORK 


A. Conclusions 



Fig. 4: Upper bounds on induced £ 2 -norms for perturbation 
velocities of Couette flow for different Reynolds numbers. 


class of appropriate Lyapunov/storage functionals for such 
flows. Conditions based on matrix inequalities are given for 
streamwise constant flows. When the laminar flow is given by 
a polynomial of spatial coordinates, the matrix inequalities 
can be checked using convex optimization. For illustration 
purposes, we applied the proposed method to study a model 
of rotating Couette flow. 

B. Future Work 

In this study, we considered flows in the Cartesian coor¬ 
dinate system. For many flows, like pipe Poiseuille flow, the 
coordinate system is naturally cylindrical. An extension of 
the results proposed in this paper to cylindrical coordinates 
is under study. 

In several scenarios in fluid mechanics, we are interested 
in a functional of the perturbation dynamics. For example, in 
the drag estimation problem, we are interested in estimating 
the functional of pressure over the surface of an airfoil. We 
are currently applying the methodology proposed in [23] to 
address such problems in fluid mechanics. 

VII. ACKNOWLEDGMENTS 


We studied stability and input-state properties of fluid 
flows with invariance in one direction. We formulated a 


The authors appreciate stimulating discussions by Prof. 
Sergei Chernyshenko from Department of Aeronautics, Im- 



































Fig. 5: Estimated critical Reynolds numbers Rec and ISS 
Reynolds numbers Reiss in terms of Ro. 


perial College London and Prof. Charles Doering, from 

Department of Mathematics, University of Michigan. 

References 

[1] P. G. Drazin and W. H. Reid, Hydrodynamic Stability. New York: 
Cambridge University Press, 1981. 

[2] D. D. Joseph, Stability of fluid motions, ser. Springer Tracts in Natural 
Philosophy. Berlin: Springer-Verlag, 1976. 

[3] J. Serrin, “On the stability of viscous fluid motions,” Arch. Ration. 
Mech. Anal., vol. 3, pp. 1-13, 1959. 

[4] V. Romanov, “Stability of plane-parallel Couette flow,” Functional 
Analysis and Its Applications, vol. 7, no. 2, pp. 137-146, 1973. 

[5] N. Tillmark and P. H. Alfredsson, “Experiments on transition in plane 
Couette flow,” Journal of Fluid Mechanics, vol. 235, pp. 89-102, 1992. 

[6] L. N. Trefethen, A. E. Trefethen, S. C. Reddy, and T. A. Driscoll, 
“Hydrodynamic stability without eigenvalues,” Science, vol. 261, no. 
5121, pp. 578-584, 1993. 

[7] P. J. Schmid, “Nonmodal stability theory,” Annual Review of Fluid 
Mechanics, vol. 39, no. 1, pp. 129-162, 2007. 

[8] P. J. Goulart and S. Chernyshenko, “Global stability analysis of fluid 
flows using sum-of-squares,” Physica D: Nonlinear Phenomena, vol. 
241, no. 6, pp. 692 - 704, 2012. 

[9] S. Chernyshenko, P. Goulart, D. Huang, and A. Papachristodoulou, 
“Polynomial sum of squares in fluid dynamics: a review with a 
look ahead,” Royal Society of London. Philosophical Transactions A. 
Mathematical, Physical and Engineering Sciences, vol. 372, no. 2020, 
2014. 

[10] O. Reynolds, “An experimental investigation of the circumstances 
which determine whether the motion of water shall be direct or sinuous 
and the law of resistance in parallel channels,” Philos. Trans., vol. 935, 
no. 51, 1883. 

[11] B. F. Farrell and P. J. Ioannou, “Stochastic forcing of the linearized 
Navier-Stokes equations,” Physics of Fluids A, vol. 5, no. 11, pp. 2600- 
2609, 1993. 

[12] B. Bamieh and M. Dahleh, “Energy amplification in channel flows 
with stochastic excitation,” Physics of Fluids, vol. 13, no. 11, pp. 
3258-3269, 2001. 

[13] M. R. Jovanovic and B. Bamieh, “Componentwise energy amplifi¬ 
cation in channel flows,” Journal of Fluid Mechanics, vol. 534, pp. 
145-183, 2005. 

[14] M. R. Jovanovic, “Modeling, analysis, and control of spatially dis¬ 
tributed systems,” Ph.D. dissertation, University of California, Santa 
Barbara, 2004. 

[15] D. F. Gayme, B. J. McKeon, B. Bamieh, A. Papachristodoulou, 
and J. C. Doyle, “Amplification and nonlinear mechanisms in plane 
Couette flow,” Physics of Fluids, vol. 23, no. 6, 2011. 

[16] M. Jovanovic and B. Bamieh, “The spatio-temporal impulse response 
of the linearized Navier-Stokes equations,” in Proceedings of 2001 
American Control Conference, vol. 3, June 2001, pp. 1948-1953. 


[17] G. Valmorbida, M. Ahmadi, and A. Papachristodoulou, “Semi-definite 
programming and functional inequalities for distributed parameter 
systems,” in 53rd Conference on Decision and Control, Los Angeles, 
CA, 2014. 

[18] M. Ahmadi, G. Valmorbida, and A. Papachristodoulou, “Input-output 
analysis of distributed parameter systems using convex optimization,” 
in Decision and Control (CDC), 2014 IEEE 53rd Annual Conference 
on, Dec 2014, pp. 4310^315. 

[19] D. D. Joseph and W. Hung, “Contributions to the nonlinear theory 
of stability of viscous flow in pipes and between rotating cylinders,” 
Archive for Rational Mechanics and Analysis, vol. 44, no. 1, pp. 1-22, 
1971. 

[20] L. Payne and H. Weinberger, “An optimal poincare inequality for 
convex domains,” Archive for Rational Mechanics and Analysis, vol. 5, 
no. 1, pp. 286-292, 1960. 

[21] J. B. Lasserre, Moments, Positive Polynomials and Their Applications. 
Imperial College Press, London, 2009. 

[22] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge 
University Press, 2004. 

[23] M. Ahmadi, G. Valmorbida, and A. Papachristodoulou, “Barrier func¬ 
tionals for output functional estimation of PDEs,” in 2015 American 
Control Conference, July 2015. 

[24] P. Parrilo, “Structured semidefinite programs and semialgebraic ge¬ 
ometry methods in robustness and optimization,” Ph.D. dissertation, 
California Institute of Technology, 2000. 

[25] M. Choi, T. Lam, and B. Reznick, “Sums of squares of real polyno¬ 
mials,” in Symposia in Pure Mathematics, vol. 58, no. 2, 1995, pp. 
103-126. 

[26] G. Chesi, A. Tesi, A. Vicino, and R. Genesio, “On convexification 
of some minimum distance problems,” in 5th European Control 
Conference, Karlsruhe, Germany, 1999. 

[27] A. Papachristodoulou, J. Anderson, G. Valmorbida, S. Prajna, P. Seiler, 
and P. A. Parrilo, SOSTOOLS: Sum of squares optimization toolbox for 
MATLAB, http://arxiv.org/abs/1310.4716, 2013, avail¬ 
able from http : //www. eng. ox. ac. uk/control/sostools. 

Appendix 

A. Poincare Inequality 

Lemma 1 ([20]): Assume H C K 2 is a bounded, convex, 
Lipschitz domain with diameter D, and u £ C 1 (Q) with no¬ 
slip u | qci = 0 or periodic such that J j n u dU = 0 boundary 
conditions. Then, the following inequality holds 

< IIVtx|| £ 2. 

B. Sum-of-Squares Programming 

Denote the ring of polynomials with real coefficients by 
R\x\, and the ring of polynomials with a sum-of-squares 
decomposition by S[i] C TZ\x\. A polynomial p(x ) £ E[a;] 
if 3 Pi(x) £ R[x], i £ {1, such that p(x) = 

H=i Pi( x )- Hence, p(x) is clearly non-negative. The set of 
polynomials is called SOS decomposition of p(x). 

The converse does not hold in general, that is, there exist non¬ 
negative polynomials which do not have an SOS decompo¬ 
sition [24]. To test whether an SOS decomposition exists for 
a given polynomial, one can solve an SDP (see [25], [24], 
[26]). SOSTOOLS [27] is a software package for solving 
SOS programs. 









